Landscape equivalent of the shoving model 
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It is shown that the shoving model expression for the average relaxation time of viscous liquids 
follows largely from a classical "landscape" estimation of barrier heights from curvature at energy 
minima. The activation energy involves both instantaneous bulk and shear moduli, but the bulk 
modulus contributes less than 8% to the temperature dependence of the activation energy. This 
reflects the fact that the physics of the two models are closely related. 
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The physics of highly viscous liquids approaching the 
calorimetricglass transition continue to attract attention 
HI HL [J H> IE HI A major mystery surrounding 
these liquids is their non-Arrhenius behavior. If r is the 
average relaxation time and r the microscopic time (of 
order 10 _13 s), the temperature-dependent activation en- 
ergy AE(T) is defined [13,111113 b Y 

Only few viscous liquids show Arrhenius temperature de- 
pendence of the average relaxation time, i.e., have con- 
stant AE(T). Most liquids have an activation energy 
which increases upon cooling. 

The "shoving" model [13J starts from the standard 
picture of a viscous liquid: At high viscosity almost 
all molecular motion goes into vibrations around poten- 
tial energy minima. Only rarely do rearrangements take 
place which move molecules from one to another mini- 
mum. This view was formulated already by Kauzmann 
in his famous 1948-review 0], and it was the starting 
point of Goldstein's "potential energy picture" Il5| which 
was recently confirmed by computer simulations jl6| . On 
the short time scale of the barrier transition - expected 
to last just a few picoseconds - the surrounding liquid be- 
haves as a solid with bulk and shear moduli equal to the 
instantaneous (i.e., high frequency) bulk and shear mod- 
uli. Just as in free volume theories the shoving model 
assumes that molecular rearrangements only take place 
when a thermal fluctuation leads to extra space being cre- 
ated locally. One may think of the surroundings as being 
shoved aside, although this cause-effect reasoning violates 
time-reversal symmetry. If there is spherical symmetry, 
the surroun ding s are subject to a pure shear displace- 
ment |7J, ll3L Il7| . Thus the work done is proportional to 
the instantaneous shear modulus Goo and one finds 
(where V c is by assumption temperature independent) 

AE(T) = Goo(T) Vc (2) 
This expression 01 fits data for the non-Arrhenius be- 



havior of several molecular liquids 13] ; in combina- 
tion with the Tool-Narayanaswami formalism the shov- 
ing model has been applied also to structural (i.e., non- 
linear) relaxations [19j. It is the short-time shear modu- 
lus which appears in Eq. because the transition itself 
is fast (compare, e.g., the analogous appearance of short- 
time friction constant in the Grote-Hynes theory [20j]). 

The basic assumptions of the shoving model may be 
summarized as follows [7J: 

• The main contribution to the activation energy is 
elastic energy. 

• The elastic energy is located in the surroundings of 
the rearranging molecules. 

• The elastic energy is shear energy. 

The purpose of this note is to give an alternative justi- 
fication of Eq. (J2J) and discuss the interrelation between 
the two approaches. 

First, we review a classical argument estimating the 
height of the barrier between two potential energy min- 
ima from the curvature around the minima. This ar- 
gument is used, e.i* in the Marcus theory for electron 
transfer reactions |2l|,|22j; in 1987 Hall and Wolynes ap- 
plied this reasonin g in their theory of free energy barriers 
in glasses 0, US • Consider first a one-dimensional 
situation with two minima the distance 2a apart (Fig. 
1). The thin curve is the potential estimated by second 
order expansions around the minima. In this approxima- 
tion the barrier height is determined by the curvature of 
the potential at the minimum: If the potential here is 
U{x) = (A/2)x 2 , the barrier is given by AE = (A/2)a 2 . 
Statistical mechanics implies (x 2 ) — fc^T/A. Therefore, 
comparing different situations with differing potentials 
but same distance between minima, the vibrational mean- 
square ir-nuctuation determines the barrier as follows: 

AE(T) « ( 3 ) 

{x 2 ) 

The basic assumptions are that the minima are at fixed 
distance and that the parabola approximation gives a 
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Configuration space 

FIG. 1: Two potential energy minima. The full curve is the 
potential, the thin curve is the potential estimated from the 
curvature at the minima. According to a classical reasoning 
used, e.g., in the Marcus theory for electron transfer reactions, 
the barrier height may be estimated by the parabolas obtained 
by extrapolating from the minima by second order expansions. 
For simplicity we show the case of identical energy minima; 
the reasoning applies, however, equally for the case of differing 
minima. 

good estimate of the barrier height. Actually, the con- 
stant of proportionality in Eq. is of no importance 
and the true barrier may be consistently smaller than es- 
timated - Eq. @ applies as long as the barrier scales 
with that estimated by the parabola approximation. 

How does the above reasoning apply to the multidi- 
mensional configuration space where the hills and valleys 
of the energy landscape live? At low temperatures a vis- 
cous liquid thermally "populates" only deep minima in 
configuration space. A transition between two deep min- 
ima most likely consists of a whole sequence of transitions 
between intermediate shallow minima Nevertheless, 
there is one or more bottlenecks in this sequence. Our ba- 
sic assumption is that at different temperatures the bot- 
tleneck transitions are the same type of local rearrange- 
ments and thus with (virtually) same distance between 
the two minima in configuration space. They occur, how- 
ever, as temperature changes in different surroundings, so 
the minima involved are different and have differing (x 2 ). 
The relevant mean-square cc-fluctuation is to be taken in 
the direction between the minima. The final assumption 
needed is that this quantity is typical for the minima, 
i.e., it is equal to the average over all directions. Con- 
sequently, the relevant mean-square x-fluctuation may 
be evaluated simply as the vibrational mean-square fluc- 
tuation averaged over all directions. We estimate this 
quantity by its ensemble average over all minima, (x 2 ). 
Note that, since the configuration space average (a; 2 ) is 
equal to the single atom mean square displacement in 
three dimensions, this quantity is experimentally acces- 



sible in, e.g., neutron scattering where it enters into the 
Debye-Waller factor exp(— Q 2 (x 2 )) (where Q is scatter- 
ing vector) [2^|. - Recently, Starr, Sastry, Douglas, and 
Glotzer arrived at Eq. from the free volume perspec- 
tive |27|. Their work includes a direct numerical confir- 
mation for a glass-forming polymer melt by calculating 
the free volume Vf and subsequently showing, with Vf as 
"mediator," that the temperature dependencies of (x 2 ) 
and r are consistent with Eq. 

The next step is to relate (a; 2 ) to the instantaneous 
bulk and shear moduli. On short time scales a viscous 
liquid behaves like a solid 0, 0, 0, |2j| . In particular, 
it has well-defined vibrational eigenstates. We assume 
that the entire phonon spectrum scales with the long 
wavelength limit of the phonon dispersion relation. For a 
one-dimensional solid there is only one elastic constant, 
C. The above A is proportional to C and consequently 
(x 2 ) oc T/C. In a three dimensional isotropic solid there 
are two elastic constants. For each wave vector there are 
three phonon degrees of freedom, two transverse and one 
longitudinal. The relevant transverse elastic constant is 
the isothermal {— adiabatic) shear modulus, while the 
relevant longitudinal elastic constant is the isothermal 
adiabatic) longitudinal modulus M defined [2^| by 
M = K + (4/3)G where K and G are bulk and shear mod- 
uli. Averaging over the two types of phonons one finds 
that the vibrational mean-square fluctuation is given by 
(x 2 ) oc T(2/G+l/M). This applies for a solid. It ap- 
plies at short times for a viscous liquid as well, if G and 
M are identified with the liquid's instantaneous isother- 
mal moduli. Thus inserting the expression for (x 2 ) into 
Eq. © we get 



AE(T) Goo(T) Moo(T) ' v ' 

Unfortunately, it is not possible to test Eq. directly 
because there are no measurements of the instantaneous 
isothermal bulk modulus. 

There are several ways to quantify variations in 
AE(T), a liquid's "fragility" [H- The stan- 
dard approach utilizes the quantity to introduced 
by Plazek, Ngai, Bohmer, and Angell [3]]]: m = 
dlog 10 (r)/(i(r s /T)| T=Tg where T g is the calorimetric 
glass transition temperature defined by r(T g ) — 10 3 s. 
Simple Arrhenius behavior corresponds to to = 16; most 
glass-forming liquids have fragilities between 50 and 150. 
As an alternative Tarjus, Kivelson, and coworkers pro- 
posed to measure the degree of non- Arrhenius behav- 
ior at any given temperature by the normalized activa- 
tion energy: AE{T) / AElT^oo) [13 It is not obvious 
a priori, however, that scaling to the high-temperature 
limit is physically relevant. Inspired by the Griineisen 
parameter (2|| and Granato's recent work on intersti- 
tialcy relations for the viscosity we suggest that a 
useful unbiased measure of how much activation energy 
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changes with temperature is its logarithmic derivative, 
din AE(T)/d\sxT, Because activation energy increases 
as T decreases, it is convenient to change sign. We thus 
define the temperature index Jab of the activation energy 
by 



dlnAE(T) 
d\nT 



(5) 



It is straightforward to show that m = 16 (l + IabC^c,))- 
We proceed to express 7a_e in terms of temperature 
indices of instantaneous moduli (same definition). First, 
note that the temperature index / of a sum, fx + / 2 , is a 
convex combination of the temperature indices h of jy. 
I = atil\ + a 2 l2 where otj = fj/(fi + fz)- Equation Q 
thus implies 



M _1 



2G^ + M^ 1 ' 2G^ + M^o 1 



Similarly, = + (4/3)Goo implies 

Im ~ = K - + 3M^ /G ° 

Substituting Eq. Q) into Eq. © leads to 



(6) 



(7) 



Iai 



(1 - a) 7 Goo + a I 



(8) 



where 



We conclude that simplifying [3_4J, but not unreason- 
able, assumptions in the landscape approach lead to a 
prediction for the non-Arrhenius behavior which in prac- 
tice is going to be hard to distinguish from that of the 
shoving model. The really interesting question is: What 
is the relation between the physics of the two approaches? 
At first sight they seem quite different. There are, how- 
ever, similarities leading us to conclude that the physics 
are actually closely related: The first shoving model as- 
sumption (main contribution to the activation energy is 
elastic energy) is equivalent to the landscape assumption 
of a parabolic potential with a curvature which deter- 
mines the activation energy. The second shoving model 
assumption (elastic energy is in the surroundings of the 
reorienting molecules) is consistent with the landscape 
assumption that the vibrational mean-square fluctuation 
is typical, because vibrational eigenstates in an isotropic 
solid involve all atoms. And the final shoving model as- 
sumption (elastic energy is mainly shear energy) is con- 
sistent with Eq. (|10|l . - Finally, we would like to re- 
mind that the original shoving model derivation of Eq. 
P|) assumed spherical symmetry |l3| . In more realis- 
tic scenarios there must be some volume change in the 
surroundings of the reorienting molecules and thus some 
contribution to the activation energy from the instanta- 
neous bulk modulus. While it is difficult to give absolute 
bounds on the magnitude of the bulk contribution, it is 
noteworth y th at Granato - supported by many others 
(see Ref. |35| and its references) - finds that for defect 
creation in a crystal the work is overwhelmingly that of 
a shear deformation. 
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2G^ + M^-, 

It is straightforward to show |3^| that one always has 



(9) 



a < 0.08. 



(10) 



Thus more than 92% of the temperature index of the acti- 
vation energy derives from the instantaneous shear mod- 
ulus. The minute influence of the bulk modulus comes 
about because of three factors: 

• There are two transverse phonon degrees of free- 
dom, but only one longitudinal. 

• Longitudinal phonons are associated with a larger 
elastic constant than transverse phonons and thus 
give less than one third contribution to the activa- 
tion energy (Eq. |@J). 

• Gqo affects also the longitudinal phonons. 
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